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' Abstract 

\ Field equations with time and coordinates derivatives of noninteger order are derived 

o 
l> 
o 

::a 

^1^' of the Ginzburg-Landau and nonhnear Schrodinger equations. As another example, dy- 

\ namical equations for particles chain with power-law interaction and memory are consid- 

ered in the continuous limit. The obtained fractional equations can be applied to complex 
media with/without random parameters or processes. 



from stationary action principle for the cases of power-law memory function and long- 
range interaction in systems. The method is applied to obtain a fractional generalization 



1 Introduction 

From the contemporary vision of complex media, where microscopic processes take place and 
where many important applications are utilized, it is too far from considering the media as 
uniform gases, liquids, or solids. The most typical features of the new physical objects and/or 
processes are fractality of their structure and of intrinsic dynamics or kinetics. Observation 
of fractality of the basic processes began fairly long ago (see for review [H [2])- Typically 
the complexity of systems is linked to the long term memory, long-range interactions, non- 
markovianity of the kinetics, and particularly with the Levy-type processes (Levy flights) [5]. 
The literature on this subject is vast. Let us mention some of the most related references. 
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where the indication of the complexity can lead, in one or another way, to the fractional 
description of the dynamic and/or kinetic processes with fractional time [H El [6]; systems 
of many coupled elements [H [8]; colloidal aggregates and chemical reaction medium (Qj [TO] : 
wave processes [HI [121 [13] ; porous media [H]; quantum mechanics and quantum field theory 
[T5l [T6l IT7] : plasma physics [HI [191 120]; magnetosphere [21]; random processes and random 
walks [221 1231 [211 [23 [26] ; fractional diffusion and Brownian motion [211 [23 [6] ; weak and strong 
turbulence [121 [23 [30] ; fractional kinetics and chaos theory [31j (see for review [321 [33] ). 

It seems that the basic formal tool to be applied to is the description of the processes by 
fractional equations, i.e., by the ones that contain fractional derivatives or integrals [3 ^ [35 1 l36 l 
[37] . The theory of derivatives of non-integer order goes back to Leibnitz, Liouville, Riemann, 
Grunwald, and Letnikov [341 [35] . Derivatives and integrals of fractional order have found 
many applications in recent studies in physics because of their continually growing numerous 
applications. 

Usually, onset of fractional derivatives (integral) is linked to different power type asymptotic 
interactions or time memories. Depending on what kind of specific features characterize the 
physical object, the fractional derivative (integral) can be with respect to time or space coor- 
dinate. In the description of particles transport, when the dynamics is chaotic, the fractional 
derivatives emerge in space and time simultaneously as a natural reflection of scaling proper- 
ties of the phase space dynamics [3ll [32] . The diffusion described by the fractional equations 
is called the anomalous one. The occurrence of such derivatives could also be related to the 
space-time decay [381 [39] , i.e., to pure dynamical processes without kinetics or diffusion. Par- 
ticularly it was shown in [13 [lH [121 23 E] how the long-range interaction between different 
oscillators can be described by the fractional differential equations in the continuous medium 
limit. Another way to connect the fractional equations with specific dispersion laws of the 
media was considered in [T21 HH US] . 

The goal of this paper is to provide a systematic approach to the onset of fractional equations 
as a result of existence of long-range interaction in a corresponding space and long-range time 
memory in the system of fields or particles depending on what kind of physical objects are 
considered. The notions of long terms memory or interaction can be exactly specified by power 
laws in time for a memory function and power law interaction between different elements of 
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the medium. It is of importance to understand the conditions when the fractional derivatives 
(integrals) occur since it allows us to involve into the consideration power tools of fractional 
calculus. 

In Sec. 2, we consider the variation of action functional that describes field with memory and 
long-range interaction. The long-time memory and long-range interaction can be introduced 
through power-like kernels of the action functional. The corresponding powers are defined 
by the exponent a (for space) and /3 (for time), which in general can be fractional. The 
Euler-Lagrange equations lead to the equation with fractional (a, /9)-derivativcs. In Sec. 3, 
the obtained results are used for derivation of (a, /3)-generalization of the Ginzburg-Landau 
and nonlinear Schrodinger equations. In Sec. 4, we consider chains of particles with long-range 
interaction and memory function. Applying the results of Sec. 2, we derive the continuous limit 
of the particle dynamics equations. In two Appendices, we provide a brief information on the 
Riemann-Liouville, Caputo and Riesz fractional derivatives used in paper, and n-dimensional 
generalization of the final fractional equations. 



2 Action functional and its variation 

2.1 Action functional 

Let us define the action functional as 

S[u] = J (fx j (fy (^dtu{x)go{x,y)dt'u{y) + ^dru{x)gi{x,y)dr'u{y) - V {u{x) , u{y))^ . (1) 

Here x — {t,r), t is time, r is coordinate, and y — {t',r'). The integration is carried out over 
a region R of the 2-dimensional space to which x belong. The field u{x) is defined in a 
2-dimensional region R of M^. We assume that u{x) has partial derivatives 

which are smooth functions with respect to time and coordinate. 
Here are three examples of this action, 
(a) If 

9o{x,y) = -gi{x,y) = S{x-y), 



V{u{x),u{y)) = V{u{xmx - y), (2) 

then we get the usual action 

S\U] = l^d'x Q [dtU{x f - \ \drU{x)f - V{U{X))^ . 

(b) If 

9o{.x,y) = -gi{x,y) = 5{x - y)Ci{D,r), 
V{u{x),u{y)) = V{u{xmx - y)Ci{D, r), 

where 

C^i(Ar) = ^L_, (0<D<1), 

then we obtain 

S[u] = jdtj dh Q [OM^)? - \ [dMx)f - V{u{x)) 
where 

dlo = Ci{D,r)dr. 

This action defines the field u{x) in a medium with the fractional Hausdorff dimension D |49] . 

(c) If 

V{u{x),u{y)) = V{u{x))5{x - y), 
9oix, y) = go6{r - r')/Co(t, t'), 

giix,y) = gi6it-t')ICi{r,r'), (3) 

then it follows from ([1]) and ([3]), 

S[u] = -go I dr I dt I dt' dtu{t,r)}Co{t,t')dt'u{t' ,r) + 
2 JR Jr Jr 

+ -91 dt dr dr' drU{t,r)}Ci{r,r')dr'u{t,r') - dt drV{u{t,r)), (4) 
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and the time and space dependent kernels are separated in the terms with derivatives. 
We will be interested in a homogeneous case 

/Ci(r,r') = /Ci(r - r'), 



and an algebraically decaying kernel /Ci with a power tail, i.e., 

/Ci(Ar) = (A)^-"/Ci(r), (1< a < 2). (5) 

Similarly, we can consider 

lCo{t,t') = ]Co{t-t') 
for < t' < t as a homogeneous function of order 1 — f3: 

/Co(At') = A^-^/Co(t'), (0 < /3 < 2, 0<t' <t). (6) 

The relation ([5]) means that we have power-law long-range interaction in the system. Equation 
([6]) indicates the memory effects with power-law memory function, which can be regarded as 
the influence of the environment. Just this case of the power-law dependences of /Co(^) and 
([5]) and IQ, will be considered to derive the field equations with fractional derivatives. 

2.2 Gateaux differential and variation of action 

The field equations will be derived by using the Gateaux differential [IHl US, HH] of S[u] at the 
point u{x), which is defined as the limit 

6S[u, h] = ( ^S[u + eh]) = hm ^l^ + e h] - S[u] ^ 

and which exists for fairly smooth integrable functions h{x) = 6u{x). The Gateaux derivative 
is slightly different from the Frechet derivative 6fS[u, h], where 

jjj^ II S[u + h]- S[u] - 6fS[u, h] \\ ^ ^ 

\\h\\^0 \\h\\ 

The Gateaux derivative is more general concept than Frechet derivative. If a function is Frechet 
differentiable, it is also Gateaux differentiable, and 6S[u,h] is a linear operator. However, not 
every Gateaux differentiable function is Frechet differentiable. In general, unlike other forms 
of derivatives, the Gateaux derivative is not linear with respect to h{x). 
The action ([1]) for u + eh is 

S[u + eh] = j d^x j d^y(^^dt{u{x) + eh{x))go{x,y)dt,{u{y) + eh{y)) + 



+^dr{u{x) + eh{x))gi{x,y)dr'{u{y) + eh{y)) - V{u{x) + eh{x),u{y) + eh{y))y 
This expression up to the order e has the form 

S[u + eh] ^ S[u] + e J (fx J fy(^^dth{x)go{x,y)dt'u{y) + ^dtu{x)go{x,y)dt'h{y)+ 

+^drh{x)gi{x,y)dr'u{y) + ^dru{x)gi{x,y)dr'h{y)- 

dV{u{x),u{y)) ^^^^ dV{u{x),u{y))^^^ .^\ ^ 



1 av[u[x ,u[y \ 



du{x) du{y) 
In the second, fourth and 6th terms of the right hand side, we change the variables x <-> 
Then 

S[u + eh] = S[u] + s I (fx J (fy(^^dth{x)[go{x,y) + go{y,x)]dt'u{y)+ 



+ ^drh{x)[gi{x, y) + gi{y, x)]dr'u{y)- 
d[V{u{x),uiy)) + V{u{y),u{x))] 



du{x) 

It is convenient to introduce the functions 



h{x)^ + ... 



Ko{x, y) = ^[go{x, y) + go{y, x)], 

Ki{x, y) = ]^[gi{x, y) + gi{y, x)], 
U{u{x),u{y)) = V{u{x),u{y)) + V{u{y),u{x)). 
Then the variation of action is 

6S[u, h] = hm ^^"^ + ^ ^1 ~ '^N ^ / f fy(dMx)Ko(x,y)dMy) + 

£ Jr Jr ^ 

+drh{x)Ki{x,y)dr'u{y) h{x)j. 

Using the relations 

dth(x)Ko(x, y)dt'u(y) = dt [h(x)Ko(x, y)dt'u(y)] - dt [Ko(x, y)dt'u{y)] h(x), 

drh{x)Ki{x, y)dr'u{y) = dr [h{x)Ki{x, y)dr'u{y)] - dr [Ki{x, y)dr'u{y)] h{x), 

dtdt'u{y) = drdr'u{y) = 0, 
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and the boundary condition 

[h{y)]dR = 0, 

we get 

dU{u{x),u{y)) 



5S[u,h]= / d'xh{x) / d'y [-dt[K,{x,y)]dMy) ~ dr[K^{x,y)]dMy) 



R JR 

For the symmetric potential 

U{u{x),u{y)) = U{u{x))5{x - y), 

equation ( ITOi) transforms into 
6S[u,h] = - [ d^xhix)( [ d^ydtKoix,y)dMy)+ f d^y drK,ix,y)dMy) + ^^^''^''^^ 



du{x) 

(10) 



Ft \J R J R 



du{x) 



(11) 

The dynamical equation follows from the stationary action principle 

5S[u, /i] = 

for any h. The field u = u{x), which leads to a minimum or saddle values of S[u], describes the 
space-time evolution. For the action ([1]), the stationary principle gives 

[ d'ydtKo{x,y)dMy)+ ! d'y drK,{x,y) dMy) + ^^^^^ = 0. (12) 
Jr Jr du{x) 

It is an integro-differential equation, which allows us to derive field equations for different cases 
of the kernels Ko{x,y) and Ki{x,y). 

2.3 Special cases 

Let us consider here two special cases: (a) system without memory and with local interaction 
in space, (b) field with power-law memory and long-range interaction. 

(a) In absence of memory and for local interaction the kernels (j9]) are defined at the only 
instant t and point r, i.e., 

Koix.y) = go6{x - y), Ki{x,y) = gi6{x - y) 
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with some constants go and gi. Then equation (1121) gives 

godMt, r) + g,dlu{t, r) + ^^^^f'"^^ = 0. 

For = 1, 5'i = -1, and 

U{u(t,r)) = — cosu(t,r), 

we get the sine- Gordon equation 

d^uit, r) - dluit, r) + sin u{t, r) = 0. (13) 

(b) In this example, we show how time and space variables can be separated leaving a 
possibility to consider the system with power-law memory and long-range interaction. Let 
Kq^x.u) and Ki{x,y) have the form 

Ko{x,y)=6{r-r')}Co{t,t'), (14) 

K,{x,y) = 6{t-t')ICi{r,r'), (15) 
where x = (t,r), and y = (t',r'). Then field equation (fT2!) can be presented as 

Z,(t,r) + Z.(t,r) + ^^M^ = 0, (16) 

where 

/ + CO 
dt'dt)Co{t,t')dMt',r), (17) 
-oo 

Zr{t,r) = dr'dr]C^{r,r')^^^^ (18) 

with separated spatial and temporal kernels. Till now, the kernels ]Co{t,t') and /Ci(r, r') were 
not defined. Their specific choice to present a long-term memory and long-range interaction 
will be in the next two subsections. 

2.4 Power- law memory 

Consider the kernel dt]Co{t,t') of integral (fT71) in the interval t' G (0,t) such that 

M(t-t'), 0<t'<t; 
dMt,t') = { (19) 

0, t' >t, t' < 0. 
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Then 



Zt{t,r)= / dt' M{t-t')dt'u{t\r) = M{t)*dtu{t,r). (20) 





As the result, we have the evolution field equation in which the quantity Zt{t,r) is related 
to another quantity dt'u{t',r) through a memory function Ai{t). Equation fl2U]) is a typical 
non-Mar kovian equation obtained in studying of systems coupled to an environment, where 
environmental degrees of freedom being averaged. For a system without memory, we have 

M{t-t') =6{t-t'), (21) 

and 

Ztit,r)= [ 5{t-t')dtm{t\r)dt' = dtu{t,r), (22) 
Jo 

i.e., the function Zt{t,r) is defined by dtu{t,r) at the only current instant t. 
Consider now the power-like memory function 

^('-0- r(i-^)(i-\r - 

where Qq is a constant that can be presented as a strength of perturbation induced by the 
environment, and r(l — (3) is the Gamma function. 
Substitution of ([23]) into ([20]) gives 

Zt{t,r) = f^^^ I (t - tr^dMt\r)dt' = go ^D^u{t,r), {0 < (3 < 1), (24) 

where f Z)f is the left fractional Caputo derivative [36l [37] . 
For the kernel dt}Co{t,t') in the integral (ITTll such that 

{M'(t'-t), t<t'<0; 
(25) 
0, t' > 0, t' < t, 

where 
we get 

^^(t,r) = ^^J^'^0^dt' = g',?D^,u, (0 < /? < 1), (27) 
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which is the right fractional Caputo derivative [371 136] . 

In general, the kernel K,Q{t,t') can include positive and negative intervals of time. Then 



dtK,^{t,t') = { 



M{t-t'), 0<t'<t; 
M'it'-t), t<t'<0; 
0, 0<t<t', t' < t < 0, 



:28i 



where Ai{t~t') and M.'(t' — t) are defined by (|23|1 and fl26l) . Then, we get a linear combination 
of left and right Caputo derivatives 



Zt{t,r) = go ^'D^u{t,r)+g', '^D^,u{t,r), (0 < /3 < 1). 



(29) 

As a result, field equation (|T2l) consists of fractional time derivatives, and it will be written in 
Sec. 2.6. 

We also will be interested in the case when 



/Co(t,t') 



M{t~t'), 0<t'<t; 

0, t' > t, t' < 0, 



(30) 



or 



Mit-t'), 0<t'<t; 

JCo{t,t') = lM'{t' -t), t<t'<0; (31) 

0, <t<t', t' <t<0, 

(compare to ( fT9l) and ( l25l) ). with the functions TW, 7W as in (!23l) and (!26|) . Substitution of 
(fT9|) and (!25|) into (fT7|) . and integration by parts gives, similarly to and (129|) . 



Zi(t,r)=^?o?A^+'«(i,^), (0</3<l), 



(32) 



or 



Zi(t,r) = o'^A^'''w(t,r) +^7^ fD^,^\it,r), (0 < /3 < 1) (33) 

with the same field equation ( fT6l) . Depending on different kernels ( !T9l) . ( l25l) and ( l30l) . ( l3Tl) . we 
obtain field equations with different order of time derivatives (see in Sec. 2.6). 
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The Caputo fractional derivatives can be linked to fractional powers of variable s for the cor- 
responding Laplace-transformed equation. It is known [501 l5T| [36], that the Laplace transform 
of the Caputo fractional derivative is 

m— 1 

= s^v{s, r) - V s^-5-1m('')(0, r), (34) 



poo 




Jo 





s=0 



where m — 1 < (3 < m, 



d'^u{t, r) 
dti 



and f (s, r) is the Laplace transform of u{t^ r) 



(s,r) = / e~''u{t,r)dt. 



(35) 



Note that formula ( IMI) involves the initial conditions u^'^\0,r) as integer derivatives u^''\t,r) 
with respect to time. Therefore we can put the initial conditions in a usual way. The functions 
u{t, r) satisfy the condition 

/■oo 

(36) 



e \u{t,r)\ < oo. 



For < /3 < 1, Eq. has the form 



Inversion of ( !37l) gives 



poo 




Jo 





^Dfuit, r) = /" e^* [sf^vis, r) - s^-'u{0, r)] ds, 



where Br denotes the Bromwich contour. 

The final equation of u{t,r) will be written in Sec. 2.6. 



(37) 



(3^ 



2.5 Non-local interaction 

Consider the kernel /Ci (r,r') of the integral ( fTSi) as 
/Ci(r, r') = C(|r — r'|) 



1 



cos(7ra/2)r(2 - a) |r - r'|"-i ' 
that describes the power law interaction. Then, we obtain 



Zr(t,r) 



1 < a< 2) 



dr'C{\r - r'\) ' = gi^^u{t, r) 



'd\r\ 



(39) 



(40) 
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where the fractional Riesz derivative with respect to coordinates is introduced [SI [37] (see also 
Appendix 1). 

It is known [M] the connection between the Riesz fractional derivative and its Fourier 
transform 

7^ — ^-i^r, (41) 

where is defined by 

/ + 00 
/(r) e-^'^dr, (42) 
-oo 

and is an inverse Fourier transform 

1 /■+00 

/(r) = (^-V~)(r) = —J Rk) e'^^dk. (43) 

The fractional Riesz derivatives describes properties of fractal media or complex media with 
fractional dispersion law (see for example in |45j ) . 

2.6 Field equations with fractional derivatives 

Substitution of ( l29ll and (l39ll into (|T6ll gives the fractional field equation 

^o?A^n(t,r) + (?^fDgn(t,r) + (?i^^(t,r)+ ^^^^ =0, (1< « < 2, < /3 < 1). (44) 

This equation describes the field of the system with power-law memory and long-range inter- 
action. Depending on the situation, qq or g'^ could be zero or not. For example, the potential 

U{u{t,r)) = \,\t,r)+^-u\t,r) 

in Eq. fj44|) gives the fractional time-dependent generalization of the Ginzburg-Landau equation 

go^D^Mt.r)+g'^^D^^u{t,r)+g^^^u{t,r) + au{t,r)+hu\t,r) =0, (1< a < 2, < /3 < 1). 

(45) 

In the case of the time kernel fl3T|) . Eq. (jUj) is replaced by 



^0 '',Dt'u{t.r)+g', ^D^,^\{t,r) +^i^«(t, r) + = 0, (1< a < 2, < /3 < 1). 

(46) 
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This equation has increased by one the order of time derivative and can be apphed to the wave 
propagation in media with fractional dispersion law. Particularly in the case when only right 
time derivative should be used (one-directional wave propagation) and the potential is 

U{u{t,r)) = — cosu{t,r), 

Eq. (j46l) gives the fractional sine- Gordon equation 

^D?^\{t,r) - -^u{t,r) + smu{t,r) = 0, {1< a < 2, < (3 < 1), (47) 

where we put Qq = 1, Qq = and gi = —1 , and which is a generalization of (fT3ll for noninteger 
derivatives with respect to time and coordinate. 

Finally, let us simplify the notation and write down Eqs. (j44|) or (H6l) as 

^) + + U\u{t, r)) = 0, (48) 

where /dt^ stays for left, right, or both Caputo derivatives (in the latter case, the constant g 
can be different for different derivatives), and 0</5<2, 0<q;<2, and U'{u) = dU jdu. Let 
us comment that the choice of the derivative jdt^ depends on the type of initial conditions 
and the processes, and other than Caputo derivative can appear. 

In Appendix 2, we present a generalization of Eqs. (43), P6|) for the ra-dimensional coordi- 
nate case. 



3 Fractional Ginzburg-Landau equation 

Since the variable x in ([1]) can be not specified, one can apply a similar technique to other 
problems, defined by the extremum of a functional with long-range interaction. As an example, 
consider a free energy functional for a model of Ginzburg-Landau equation (GLE) that consists 
of long-range interaction. 

The fractional generalization of the Ginzburg-Landau equation (FGLE) was suggested in 
Ref. [15]. This equation can be used to describe the processes in complex media [52l[53]. Some 
properties of FGLE are discussed in [H], [SU [55] . 
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It is known [56] that the stationary GLE 

g^u — au — bu^ = 

can be derived as the variational Euler-Lagrange equation 

SF[u] 



5u{r) 

for the free energy functional 



(49) 



F[u] = Fo + - / [g{duy + au^ + -u^]dr, 

^ JR ^ 



(50) 



where du = du{r)/dr, and the integration is over a region R. Here Fq is a free energy of the 
normal state, i.e. F[u] for u = 0. 

Consider the thermodynamic potential (free energy functional) F[u] for the non-equilibrium 
state of a medium with power-law non-local interaction. The generalized free energy functional 
has the form 

F[u]=Fo+ I dr I dr' J^{u{r),u{r'),du{r),du{r')), (51) 
Jr J r 

where the generalized density of free energy 

J^{u{r)^u{r-'),du{r),du{r')) = 

= ^gK:iir,r')d'^uir)d^,uir') + (^^u\r) + ^u\r)^ 5(r - r') 
has the kernel /Ci (r,r') defined as in fl5^ . The variational equation gives 



(52) 



Qa 

g-^^u{r) + au{r) + bu^{r) = 0, (1 < a < 2), (53) 

which can be called the a-FGLE. This equation can be easily generalized for the 3-dimensional 
variable r (see Appendix 2). 

In the non- stationary case, Eq. ( 149|) should be replaced by 

Mt,r) ^ 6F[u] 

dt 6uit,ry ^ ' 

(see [57]), where there is no explicit time memory effects. To put such memory into fl54l) . we 
can write 



I dt'M{t-t 
Jo 



,.du{t',r) 6F[u{t,r) 



dt' 6u{t, r) 
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(55) 



and to assume for Ai{t—t') the power law (!23l) . Then we arrive to a nonstationary generahzation 
of (a,/3)-FGLE 

-^^ = g^u{t,r) + au{t,r) + bu'{t,r), {0<(3<1, l<a<2), (56) 

where d^/dt^ is used for Caputo derivative while any other fractional derivative can be applied 
by modifying the memory kernel M.{t), and initial conditions. 

It is worthwhile to compare Eq. (l56i) to its counterpart nonlinear Schrodinger equation 
(NSE) 

.du ^ ,2 

t— = g/lu + au + o\u\ u (57) 

with u = u{t,r), and complex a, b. In the case of A_|_ instead of A, Eq. (157|) also known as 
parabolic equation for wave propagation. Generalization of fl57|) for the case of fractional space 
derivative and non-local interaction (a-NLS) was considered in [151 155] : 

diii 

i— = -g{-A)''^\ + au + b\u\\, (1 < a < 2), (58) 
where the fractional Laplacian is defined through the Fourier transform and Riesz derivatives 

m- 

T : (-A)"/2 — > (k^)"/^ (59) 
Similar to fl56|) generalized (a,/5)-NLS equation has the form 

^ = -g{-Ar^''u + aw + V (60) 
where jdt^ is now Riemann-Liouville derivative with the Fourier transform 

(61) 

and the memory function is working through the parameter /3. 

It is convenient also to interpret fl60|) through the nonlinear dispersion law by applying to 
( !60|) Fourier transform in both time and space. Then it gives with the help of (!59|) and (!6T]) 



1/3 



-^7(k2)"/=^ + a + 6|u|, (62) 



which was derived for /3 = 1 in |35l [55]. Onset of fractional time derivative in (l60l) can stop 
self-focusing of waves, steepening of the solution, developing of a singularity. These phenomena 
need a special analysis. 
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Eq. (l60l) can be easily generalized for the anisotropic case 

^ = -g^{-^^Y-l\ - ^|,(-A||)"ii/2« + au + h\u\\ (63) 
with a corresponding anisotropic dispersion equation instead of ( !62l) (see also [45l[55] for [3 = 1). 

4 Discrete system with memory and long-range interac- 
tion 

4.1 Equation for discrete chains 

In this section, we show how the obtained results of Sec. 2 can be applied to discrete systems, 
for example chains of interacting particles. 

Long-range interaction is a subject of a great interest since a long time. Thermodynamics 
of a model of classical spins with long-range interactions has been considered in [581 EH EOj 
[6H [62] . The long-range interactions have been widely studied in discrete systems of lattices 
as well as in their continuous analogues: solitons in one- dimensional lattice with the Lennard- 
Jones-type interaction [63]; kinks in the Frenkel-Kontorova model [6l|; time periodic spatially 
localized solutions (breathers) [651 ES]; energy and decay properties of discrete breathers in 
the framework of the Klein-Gordon equation |68j, and discrete nonlinear Schrodinger equations 
[19] . A remarkable property of the dynamics described by the equation with fractional space 
derivatives is that the solutions have power-like tails. Similar features were observed in the 
lattice models with power-like long-range interactions [691 EHl [661 EOl [711 US]- Long-range 
interaction can be relevant to the systems such as neuron populations [72] and Josephson 
junctions [73]. The syncronization of chaotic systems with power-law long-range interactions 
were considered in [731 [lH [71] . A model of coupled map lattices with coupling that decays in a 
power-law, was considered in [751 [761 [771 [74] . Note that the fractional power-law dependence can 
be linked to fractal properties of heterogeneous surfaces [9] and power-law decay of structure 
factor for geometry of colloids aggregates [TD]. It will be shown how long-range coupling of 
particles and memory function with power tails can reveal a new type of particle equations 
with fractional derivatives and the connection of these equations to their continuous media 
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counterpart . 

Consider an one-dimensional chain of interacting oscillators that can be described by the 
action, 

(64) 



dt / dt' ^ C{Un{t),Un{t'),Un{t),Un{t')), 
-°° n=-oo 

where m„ are displacements of the oscillators from the equilibrium and £ is a Lagrangian. If 

C{Un(t),Un{t'), Un{t), Un{t')) = C{Un(t) , Unit)) 6 {t - t'), 

then we have the chain without memory. 

Let us introduce a generalization of flM|l with the action 

'rOO p + OO . 



S[Un] 



dt 

— oo J — oo 



dt'{j: 



ICo{t,t')Un{t)Un{t') - V (Unit), Unit')) 



- J2 UiUnit),Un,it')) 



(65) 



n,m=— oo 



In the same way as in Sec. 2, let us separate the kinetic energy from the long-range interaction 
and potential parts: 



UiUnit), Umit')) = -goJai\n - m\) iUnit) - U^mit))^ 5(t - t') 
ViUnit),Unit')) = ViUnit))6it-t'). 

Note that (1^ . (1^ in (1M|) are equivalent to 

Uiunit),Umit')) = ~goJai\n - m | ) M„ (t)Mm(t) 6it - t'), 

ViUnit),Unit'))= (Viunit)) + ^~gulit)] 5(t-t'), 



(66) 
(67) 

(68) 
(69) 



where 

9 = 9oY^ Jai\m\). 

The second term in the right hand side of fl69l) removes the infinity of the interaction fl68l) in 
the continuous medium limit. The interparticle interaction J„(|n — m|) in (166!) is defined by 

1 



Jai\n - m\) 



\n — m 



a+l ' 



{a > 0). 



(70) 
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Some other examples of functions Jain) can be found in 
Using ([19]) and ([23]), for the kernel /Co(t,t') in ([65]), i.e. 

, nf^^it - tT^ ^ 0<t'<t (0</3<l); 
0, t' > t, t' < 0, 

we obtain the corresponding Euler-Lagrange equations 

+ 00 

^AV(t) +90 - "^D ["™(^) - + ^K(^)) = 0' (72) 

m=—oo 

where F{u) = dV{u)/du. 

A continuous limit of equation (1721) can be defined by a transform operation from u„(t) to 
u{x,t) ^0, HU 142] HH m]. First, define Unit) as Fourier coefficients of some function u{k,t), 
k e [-K/2,i^/2], i.e. 

+ 00 

u{t,k)= J2 «nWe-''""=-^AK(t)}, (73) 



where x„ = nAx, and Ax = 27r/i^ is a distance between nearest particles in the chain, and 

1 f+K/2 
-K/2 

Secondly, in the limit Ax {K — > oo) replace Un(t) = {2n / K)u{xn,t) u{x,t)dx, and 
Xn = nAx = 27m/ K x. In this limit, Eqs. ( l73l) . ( l74l) are transformed into the integrals 



1 r+K/2 

Unit) = - dk u{t, k) e*'^^" = J'^'iuit, k)}. (74) 

^ J-K/2 



/ + 00 
dx e-'^'^uit.x) = J^{u{t,x)} = lim J^A{un{t)}, (75) 
■oo ^^^^ 

1 /" + 00 

M(t, x) = — dk e*'=^M(t, A;) = J^-^{u{t, k)} = lim ^■^^^(t, k)}. (76) 
27r Ax^o 

Applying (1731) to (1721) and performing the limit (175|l . we obtain 

d^u{t,x) d°'u{t,x) 



dtP a|x| 

where 



+ F(M(t, x)) = 0, (0 < /3 < 2, 1 < a < 2), (77) 



^„ = 2^o(Ax)"r(-a)cos(^) {7i 
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is the renormalized constant. The Caputo time derivative is written in a simphfied form /dt^, 
and the value of (3 depends on the choice of memory function. The equation flTTl) can be 
generahzed to a nonhnear long-range interaction. Consider, instead of (172|) . 

+00 

+ 90 Y1 -^"(I^ - "^D [^(^-) - ^(^")] + ^("") = 0' ('^9) 

m,=— 00 

where /(u) is a function of m. For example, f{u) = or /(n) = u—gv^. Then the corresponding 
continuous limit for the same Ja(|n — m|) as (!70l) leads to the time-space fractional equation 

^ + 9ag^fiuit, x)) + x)) = 0, (0 < /3 < 2, 1< a < 2). (80) 

Equations ([77]) and (IHQl) for /5 = 1 were considered in (iQl El US 031 E] . Generalization to 
< /3 < 2 significantly extends the area of their applications. A physical motivation is that a 
dynamical process typically reveals fractional features simultaneously in space and time. Such 
situation just was considered in chaotic dynamics [3T1 [32] . Now we have such a possibility 
far beyond the fractional kinetics. An evident generalization of (IHOll is for the inter-particle 
interactions with two or more different kernels. For example one can consider regular terms 
without long memory together with a term with long memory: 

d^uit.x) d^u(t,x) d"u(t,x) ^/ / XX „ ^ ^ . ^x /^,x 

-^^ + 3s^^ + 9.^l^ + F{u{t,x)) = 0, (0</3<2, 1<«<2) (81) 

with some gs and ga and integer s. 



5 Conclusion 

Starting from a variation of the action functional, we consider different type of kernels that 
define the character of particle interaction and the influence of an environment on the memory 
function. The main stress is on the long-range interaction and memory that occur in complex 
media. The case when the interaction or memory function have power-law structure the system 
can be described by the equation of motion with fractional derivatives d^/dt^ and d°'/d\x\°' 
depending on the power of interaction and memory function. We have discussed how different 
types of the derivatives and possible values {a, jd) may occur with respect to the type of memory 
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and interaction. The final equations of motions can be considered as a new kind of tool to 
study dynamics with space-time distributed interactions. Number of examples of such kind of 
systems can be found in the reviews [H [32] related to random or chaotic processes. The study 
of this paper shows that the list of possible applications of fractional equations can be naturally 
expanded to include non-chaotic and non-random dynamics as well. 

This work was supported by the Office of Naval Research, Grant No. N00014-02-1-0056, 
and the NSF Grant No. DMS-0417800. 
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Appendix 1: Fractional derivatives 

The fractional derivative has different definitions [5H ES] , and exploiting any of them depends 
on the kind of the problems, initial (boundary) conditions, and the specifics of the considered 
physical processes. The classical definition is the so-called Riemann-Liouville derivative [311 ES] . 
The left and right Riemann-Liouville derivatives for an interval [a, b] are defined by 



V{n- a) dx" {x - 2)"-"+! ' 
;-l)" (9" u{z)dz 



tV?u[x) = -j / 7 -T, 82 

where n — 1 < a < n. End-points a, b can be extended to —00, 00 if the integral exists. 

Due to reasons, concerning the initial conditions, it is more convenient to use the Caputo 
fractional derivatives [36J . Its main advantage is that the initial conditions take the same form 
as for integer-order differential equations. The Caputo fractional derivatives are 

where = d'^u{z)/dz^, and n — 1 < a < n. The Caputo fractional derivatives can be 

defined through the Riemann-Liouville derivatives |37j by 



n—l , 
k=0 



C 



" ^ (6 - x)^ 



:D^u{x) = .V- I u{x) - J2 ^^«^'^(^) I ' (84) 

k=0 



where n — 1 < a < n. These equations give 



^ ^ ( \k — OL 

aVy{x) = ^D>(x) + (85) 

i (fc - a + ij 

,.P>(x) = ^^Dtu{x) + r?~''^'~], n^'\b), (86) 

^r(A;-a + l) 
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The Riesz fractional derivative of order a are 



d\xY 2cos(7ra/2) 
where a 7^ 1,3,5..., and are Riemann-Liouville fractional derivatives with infinite limits 



+ ^ ^ r(n - a) (x - ' 

" ^ ^ r(n - «) - 

Substitution of Eqs. ([HSD into Eq. dHTj) gives 

^ «w ^ f r ..^^ . r I . (89) 



2cos(7ra/2)r(n - a) (9x" Vi-oo (a^ - -s)"""'^^ (-^ - 

The Fourier transform of the fractional derivatives [Ml [37] are 



0'[^u{^)\{k) = -Wu{k). (90) 



^ (P±u(x)) (A;) = {±%kYu{}z), (91) 

where T is defined by 

/+00 
u[x) e-'^^'dx. (92) 
-00 

The inverse Fourier transform is 

«(x) = (^-1m)(x) = — / u(k)(^^^dk. (93) 

Appendix 2: n-dimensional case 

The generalization the action ([1]) for the case r G M", where x = (t, r), and r = (x^, x"), 
gives the field equation 

,0 o^A'^^t, r) + f Do^.(t, r) + |: 9.^u{t^ r) + = 0. (94) 

For the case r G M", there exists the other possibility to define the kernels Mt{x,y) and 
Mr{x, y) = {Mfc(x, y),k = 1, n}. We can consider 

Zr{t,r)= [ d-r'j2Mdr,r')^^f^, (95) 
Jr" c'x 
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where Mk{r,r') is Riesz kernel [M] : 



Mk(r,r') = KaJr - r' 



7n(afc) 



ak-n^ 0,2,4, 
r — r'|"''~" In |r — r'| — n = 0, 2, 4, 



Here > 0, (0;^ 7^ n, n + 2, n + 4, ...), and 



7n(a) = < 



2"7r"/2r(a/2)/r(^) a 7^ n + 2fc, n ^ -2A;, 

1 n = -2k, 

(_l)(«-")/22"-i7r"/2 r(a/2) [^]\ a^n + 2k. 

Note that the muhivariable Riesz integral 



where a > 0, can be presented as convolution: 

{ru) (t, r)= I Ka{r- r') u{t, r') rfV, 
with the Riesz kernel Kr,(r). It allows us to write fl95l) as 



z,(t,r) = 5^r'= 



fc=l 



r) 



(96) 



(97) 



(98) 



(99) 



The fractional Riesz integrals of orders au {k = l,...,n) in the field equations describe the 
fractal media. Then the field equation is 



k=l 



dx'^ du{t, r) 



(100) 



If Mk{r,r') in fl95p is an operator such that 



Mk{r,r')u{t,r') = gu 



1 (A^^)(t,r) 



(101) 



where 



2-a^l+n/2 



(-1)^-^/! 

r(l + a/2)r((n + a)/2) sin(a7r/2) ^ (/ - k)\k\ 
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is normalized multiplier [37], and 



n 



k=0 



is symmetrized difference [3^, then 



k=l ' ' 



As a result, we have 



which is the field equations with n fractional Riesz multivariable derivatives. 
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